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摘 要 : 准确 及 时 地 绘制 森林 地 上 碳 密 度 图 是 了 解 全 球 碳 循环 的 必要 条 件 。 虽 然 星 载 激 光 雷 达 
(如 ICESat/GLAS) 数 据 已 被 广泛 用 于 估算 大 尺度 的 森林 地 上 碳 密 度 分 布 ,但 地 形 坡 度 对 GLAS 提 取 
冠 层 高 度 精度 的 影响 往往 限制 了 其 在 山区 森林 的 应 用 。 通 过 以 祁连山 国家 公园 为 研究 区 域 ,结合 


GLAS 数据 、Landsat OLI(Operational land imagery ) 数 据 、 样 地 调查 数据 ,对 祁连山 地 区 进行 区 域 性 的 
森林 地 上 碳 密 度 估 算 。 首 先 通过 改进 后 的 地 形 校正 模型 减 小 坡度 对 GLAS 数据 提取 和 森林 冠 层 高 度 
精度 的 影响 ,使 得 更 多 的 GLAS 数 据 可 用 于 后 续 的 研究 ;其 次 将 所 建立 的 不 同类 型 森林 地 上 生物 量 
估算 模型 与 相关 碳 含 量 转 换 系 数 结合 ,得 到 GLAS 激 光 光 斑 ( 脚 印 点 ) 的 森林 地 上 碳 密 度 ; 最 后 禾 


= 


用 


FAAR He KIB (MaxEnt ) 4 2! £8 2 E L BA 2018 年 森林 地 上 碳 密度 分 布 图 。 结 果 表 
明 : 祁 连 山 国家 公园 2018 年 平均 森林 地 上 碳 密度 为 40.72+6.72 t hm?, E 1A BH 28.5824.72 Tg, 
海拔 2770~3770 m 区 域 的 森林 植被 碳 储 量 最 大 , 且 阴 坡 的 碳 储量 明显 高 于 阳 坡 。 采 用 森林 资源 清 
查 数据 独立 验证 估算 结果 的 准确 性 ,模型 估 测 均 方 根 误差 (Root mean square error, RMSE) 77 18.946 
thm?。 本 研究 结果 可 以 为 监测 区 域 乃 至 国家 尺度 的 森林 碳 储量 变化 以 及 制定 可 持续 的 森林 管理 
措施 提供 依据 。 此 外 ,本 文 所 采用 的 方法 在 山区 森林 矶 储量 估算 方面 也 具有 和 较 大 的 潜力 。 

X 键 i: 星 载 激光 雷达 ; Landsat OLI; 森林 地 上 碳 密度 ; 非 参数 化 算法 ; 祁连山 
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森林 作为 陆地 生态 系统 中 最 大 的 碳 库 ,在 维持 
全 球 碳 平衡 和 减缓 气候 变化 方面 发 挥 着 重要 的 作 
用 "。 尽 可 能 准确 地 估算 森林 地 上 碳 密度 及 其 空间 
分 布 ,是 理解 全 球 碳 循环 以 及 制定 减缓 气候 变化 相 
关 措 施 的 关键 ”。 

遥感 技术 由 于 其 具有 覆盖 范围 广 、 重 复 观 测 能 
力 强 等 特点 ,目前 已 成 为 大 尺度 森林 地 上 碳 密 度 估 
算 的 首选 。 特 别 是 激光 雷达 遥感 (Light detection 
and ranging, LIDAR) ,由 于 其 具有 直接 测量 森林 垂 
直 结 构 的 能 力 ,能 够 克服 光学 及 微波 遥感 在 高 生物 
量 地 区 出 现 的 信号 人 饱和 问题 ,已 被 IPCC(The Inter- 
governmental Panel on Climate Change ) 推 荐 为 森林 碳 
密度 估算 的 有 力 工 具 ” ,并 被 广泛 应 用 于 目前 的 森 
林地 上 碳 密度 估算 研究 “” 。 虽 然 机 载 激 光 雷 达 数 
据 能 够 提供 更 高 精度 的 森林 垂直 结构 信息 ,但 较 高 
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的 使 用 成 本 和 较 大 的 数据 处 理 量 通常 限制 了 其 在 
大 尺度 的 应 用 ""。 而 星 载 激光 雷达 如 搭载 在 冰 、 云 
和 陆地 高 程 卫星 上 的 地 球 科 学 激光 测 高 系统 
(GLAS) ,虽然 只 能 对 激光 光斑 (脚印 点 ) 照 射 范围 内 
的 植被 进行 采样 ,但 由 于 其 全 球 范围 的 采样 策略 以 
及 能 够 免费 获取 的 数据 ,目前 已 有 众多 学 者 将 其 
与 光学 遥感 数据 相 结合 对 和 森林 地 上 生物 量 / 碳 储量 
进行 估算 并 取得 了 较 好 的 结果 *"。 

采用 GLAS 数 据 估算 森林 地 上 碳 密度 的 一 个 关 
键 问题 是 如 何 准确 获取 脚印 点 内 和 森林 垂直 结构 信 
息 。 由 于 GLAS 波 形 对 地 形 较 为 敏感 , 当 坡 度 大 于 
10° 时 ,植被 高 度 的 估计 精度 明显 降低 "”。 因 此 , 许 
多 森林 碳 储 量 估 算 人 研究 去 除了 坡度 较 大 地 区 的 
GLAS 数 据 以 避免 误差 ,这 造成 了 对 山地 森林 碳 
储量 的 低估 。 而 山地 森林 由 于 其 受 人 为 干扰 相对 
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BUD ,在 维持 碳 贮存 方面 发 挥 着 关键 作用 , 且 不 同 
区 域 森 林 的 固 碳 能 力也 存在 差异 ,因此 ,提升 GLAS 
数据 在 山区 获取 的 和 森林 垂直 结构 精度 , 尽 可 能 准确 
地 估算 山地 和 森林 地 上 碳 密度 ,补充 以 往 研 究 对 山地 
森林 矶 储量 的 低估 ,对 理解 区 域 乃 至 全 球 碳 循环 具 
有 重要 的 意义 。 
祁连山 作为 我 国 西北 干旱 地 区 主要 山系 之 一 ， 

是 我 国 西部 重要 的 生态 屏障 和 固 碳 场 所 5。 目前 
针对 祁连山 地 区 森林 碳 储 量 的 估计 仅 限 于 较 小 地 
区 的 特定 树种 ””。 在 此 背景 下 ,本 研究 结合 CLAS 
数据 .Landsat OLI 数 据 和 样 地 调查 数据 ,对 祁连山 
家 公园 森林 地 上 碳 密度 进 行 估算 , 旨 在 完善 目前 
山地 区 域 的 森林 碳 储量 现状 分 析 , 为 监测 区 域 乃 至 

国家 尺度 的 森林 碳 储量 变化 以 及 制定 可 持续 的 森 
林 管 理 措施 提供 依据 。 


1 研究 区 概况 


祁连山 国家 公园 地 理 位 置 在 36*45'~39"47'N， 
94°50'~102°59'E Z [B] ,总 面积 5.02x10' km’, Hore HE 
INB 34400 km’, FYE 15800 km*( 图 1)。 该 地 区 海 
拔 1770~5740 m, 坡 度 0"~58"。 属 典型 的 温带 大 陆 
性 山地 气候 ,年 平均 降水 量 450~650 mm ,年 平均 气 
温 -3~1 CC。 区 内 森林 主要 由 青海 云 杉 (Picea crassi- 


96°E 98°E 


folia) 、 祁 连 圆 柏 (Juniperus przewalskii) 、 山 杨 Ee 
lus davidiana ) Ñl H T6 (Betula platyphylla) 组 成 。 其 
他 散 生 树 种 有 油 松 (Pinus tabuliformis) | 5€ tH BK 
(Quercus mongolica ) FU {il f£] (Platycladus orientalis ) 
等 。 区 内 灌木 树种 主要 有 高 山 绣 线 菊 (Spiraea alpi- 
na) Aib EX (Euonymus frigidus) BERZE ( Ceratoi- 
des latens) 、 金 露 梅 (Potentilla fruticosa) 等 ,草本 植物 
EZA P CES (Carex lanceolata) RAE (Polygo- 
num viviparum ) PU i LL 5 5! (Kobresia pygmaea ) 等 。 
青海 云 杉 主要 分 布 在 海拔 2500~3300 m 的 阴 坡 和 半 
明 坡 上 , 祁 连 圆 柏 主要 分 布 在 海拔 2700~3300 m 的 
阳 坡 , 半 阴 坡 或 半 阳 坡 上 ,党 与 云 杉 混交 。 山 杨 和 
a FEE UAE PR AY EO BP EZD BS 
溪流 边 和 海拔 2700 m 以 下 的 洪 泛 区 。 


2 数据 来 源 和 研究 方法 


2.1 数据 来 源 与 处 理 

2.1.1 GLAS 数据 本 研究 所 用 CLAS 数据 为 通过 
美国 国家 冰雪 数据 中 心 (https://NSIDC.org/data/ice- 
sat) 下 载 的 2 级 CLA14 产 品 。 下 载 的 详细 数据 集 
包括 L3F (2006 年 5 一 6 月 )、L3G (2006 年 10—11 
H).L3H(2007 4E 3—4 H ) .L31(2007 ^F. 10—11 H) 
L3J (2008 4E. 2—3 H ) fl L3K (2008 4E 10 HO. 研究 
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1 研究 区 位 置 示意 图 
Fig. 1 The location of the study area 
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区 内 GLAS 脚 印 点 的 空间 分 布 如 图 1 所 示 。 

下 载 后 首先 依照 以 下 2 个 规则 对 无 效 数据 进行 

过 滤 :(1) 使 用 GLA14 产品 中 记录 的 检测 标志 
(i_satNdx=0,i_FRir_qaFlag=15) 识 别 无 饱和 以 及 无 
zB GLAS BHR IFLR A?) 5 (2) 将 GLA14 数 据 中 记 
录 的 高 程 (d_elev) 与 DEM 数据 进行 比较 ,将 高 差 阐 
值 设 为 略 大 于 该 地 区 最 大 冠 层 高 度 (50 m) ,将 大 于 
此 高 差 的 数据 视 为 潜在 的 低 云 数据 并 删除 小 。 过 
滤 后 ,共有 571 个 CLAS 波形 数据 可 用 ,将 其 转换 为 
WGS84 坐标 系 以 与 其 他 数据 相 匹 配 。 
2.1.2 Landsat OLI 数据 ”本 人 研究 所 使 用 的 Land- 
sat 0LI 数 据 从 美国 地 质 调查 局 国家 地 球 资 源 观测 
与 科学 中 心 (http://GLOVIS.usgs.gov/) 获 取 , 下 载 云 
量 小 于 3% 的 Landsat LIT 数 据 共 9 景 (P131R34_ 
20180721,P132R34 20180712,P133R33 20180715, 
P133R34. 20180621, P134R33. 20180729, P135R33 
20180717,P136R32, 20180809, P136R33 20180825, 
P137R33_20180612)。 获 取 后 首先 对 其 进行 辐射 定 
标 及 大 气 校正 处 理 ,并 通过 实地 调查 建立 51 个 地 面 
控制 点 ,采用 二 次 多 项 式 数学 模型 对 影像 进行 几何 
精 校 正 , 将 误差 控制 在 0.5 个 像 元 以 内 中。 最 后 , 利 
用 直方 图 匹配 算法 对 图 像 进行 拼接 并 裁剪 以 使 其 
完整 覆盖 研究 区 域 。 
2.1.3 样 地 调查 数据 本 研究 的 样 地 调查 分 为 两 阶 
段 进行 :2018 年 10 月 及 2019 年 7 月 。GLAS 脚印 点 
的 实地 调查 于 2018 年 10 月 进行 。 由 于 不 同 地 形 坡 
度 对 GLAS 波 形 的 影响 可 能 不 同 ,因此 我 们 首先 将 
研究 区 坡度 分 为 5 级 :0°~5° 、5°~15° 、15°~25?° 、25°~ 
35°、35°~58?, 并 分 别 在 各 坡度 级 别 选取 位 于 森林 区 
域 的 GLAS 脚印 点 进行 调查 ,其 中 14 个 脚印 点 位 于 
0°~5° ,30 个 脚印 点 位 于 5°?~15° ,20 个 脚印 点 位 于 
15°~25° ,26 个 脚印 点 位 于 25°%~35° ,7 个 脚印 点 位 于 
35°~58°, 

采用 样 方 调查 法 对 GLAS 脚印 点 进行 调查 。 由 
于 本 研究 中 使 用 的 GLAS 数 据 仅 来 自 3 号 激光 器 且 
其 脚印 点 椭圆 度 最 小 并 接近 直径 为 55 m B I] 7, y 


具体 样 地 布局 如 图 2 所 示 。 使 用 精度 约 为 3 m 的 手 
持 全 球 定位 系统 (GPS) 接 收 器 确定 采样 点 中 心 坐 
标 。 用 手持 激光 高 度 计 (高 度 分 辩 率 为 0.01 m ) 测 量 
样 地 内 胸径 大 于 5 cm 的 每 棵 树 的 高 度 ,并 记录 其 胸 
径 。 同 时 记录 样 地 内 树种 、 冠 层 郁 闭 度 、 林 分 类 型 
和 地 形 。 脚 印 点 样 地 平均 冠 层 高 度 定义 为 各 子 样 
地 平均 冠 层 高 度 的 平均 值 。 


图 2 GLAS 脚印 点 样 地 布局 图 
Fig. 2 Conceptual layout of sample plot 


for measuring canopy heights in the GLAS footprints 


由 于 GLAS 数 据 的 不 连续 性 ,位 于 GLAS 脚 印 点 
内 的 森林 并 不 能 代表 祁连山 典型 林 分 。 为 建立 更 
准确 的 森林 地 上 生物 量 估算 模型 ,于 2019 年 7 月 进 
行 了 第 二 次 样 地 调查 。 此 次 调查 选取 研究 区 内 具 
有 代表 性 树种 和 林 龄 的 林 区 共 设 置 98 个 直径 为 15 
m 的 圆 形 样 地 。 记 录 样 地 内 胸径 大 于 5 em 树 的 高 
度 及 胸径 。 同 时 记录 样 地 中 心 坐标 、 树 种 、 冠 层 郁 
闭 度 、 林 分 类 型 和 地 形 。 此 外 ,在 2 次 样 地 调查 的 过 
程 中 记录 了 人 研究 区 具有 代表 性 的 景观 ,以 验证 土地 
履 被 类 型 的 分 类 精度 。2 次 样 地 调查 点 的 分 布 如 图 
3 所 示 。 
2.1.4 辅助 数据 ”本 研究 使 用 的 DEM 数据 由 研究 
区 1:50000 地 形 图 数据 生成 。 DEM 数 据 在 本 研究 中 
除了 用 来 过 滤 GLAS 数 据 以 及 消除 地 形 对 GLAS 数 
据 的 影响 之 外 ,由 其 生成 的 研究 区 海拔 .坡度 及 坡 


将 样 地 调查 样 点 分 辨 率 与 CLAS 数据 分 辨 率 相 匹 
fic ,我 们 以 所 选 CLAS 脚印 点 中 心 为 圆心 建立 直径 
为 55 m 的 圆 形 采 样 点 进行 采样 ,每 个 样 地 分 为 4 个 
圆 形 子 样 地 ,1 个 位 于 主 样 地 的 中 心 , 男 3 个 位 于 中 
心地 块 360°? .240? 和 120° 的 方位 角 , 每 个 子 样 地 的 
直径 为 15 m, 每 个 子 样 地 中 心 之 间 的 距离 为 20 m, 


向 分 布 图 还 被 用 来 作为 环境 数据 与 GLAS 数 据 相 结 
合 估算 森林 地 上 碳 密 度 的 空间 分 布 。 

本 研究 使 用 的 环境 数据 包括 从 世界 气候 数据 
网 站 (http:/www.WorldClim.ore/) 下 载 的 空间 分 辩 率 
为 1km 的 19 个 生物 气候 变量 CeoTif 文 件 ” 以 及 从 
中 科 院 资源 环境 数据 云 平 台 下 载 的 土壤 质地 空间 
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研究 区 坡度 等 级 
mum 0°~5° 
mm 5°~15° 
mum 15°~25° 
EH 259—359 
EH 359.58? 


e 2018 年 10 月 样 地 调查 点 
e 2019 年 7 月 样 地 调查 点 


96°E 98°E 


注 : 底 图 为 从 DEM 数据 中 提取 的 坡度 
图 3 样 地 ) 


40°N 


38°N 


100°E 102°E 
图 ,所 有 样 地 调查 点 均 位 于 甘肃 省 境内 。 
周 查 点 分 布 图 


Fig. 3 Distribution of field inventory plots 


分 布 数据 。 所 有 辅助 影像 数据 均 利用 3 次 卷 积 采样 
方法 将 影像 像 元 重 采 样 为 30 m。 采 用 人 研究 区 森林 
资源 二 类 调查 小 班 数 据 ( 以 下 简称 森林 资源 清查 数 
据 ) 对 森林 冠 层 高 度 估算 结果 及 森林 地 上 碳 密度 估 
算 结 果 进 行 独 立 验 证 。 

2.2 研究 方法 

2.2.1 研究 区 土地 覆盖 分 类 ”人 研究 区 的 土地 覆盖 分 
类 采用 监督 分 类 法 。 和 森林 类 型 训练 样本 来 自 于 森 
林 资 源 清查 数据 ,其 他 类 别 的 训练 样本 来 自 Google 
Earth "高 分 辨 率 影 像 , 各 类 别 训练 样本 具体 信息 如 
表 1 所 示 。 训 练 样本 确立 后 ,将 其 输入 到 最 大 似 然 
分 类 器 中 得 到 分 类 结果 ,以 实地 调查 数据 作为 验证 
样本 ,得 到 该 区 域 的 分 类 总 体 精度 为 91.59% ,Kappa 
系数 为 0.89。 通 过 分 类 结果 ,提取 位 于 森林 区 域 的 
GLAS 脚印 点 共计 273 个 ,其 中 针 叶 林 91 个 , 阔 叶 林 
55 个 , 针 冰 混交 林 127 个 。 

2.2.2 GLAS 和 脚印 点 的 森林 冠 层 高 度 估 算 在 
GLA14 产 品 中 ,记录 了 GLAS 原始 波形 通过 高 斯 滤 
波 器 拟 合 的 1 个 波形 信号 开始 范围 增 量 (d_SigBe- 
gOff) 以 及 最 多 6 个 高 斯 波峰 质心 范围 增 量 
(d. gpCntRngOff) ,本 研究 以 d_SigBegOff 作 为 植被 冠 
顶 波形 信号 ,以 最 后 2 个 d_gpCntRngOff 中 较 强 峰 近 


#1 训练 样本 详细 信息 


Tab.1 Detailed information of training samples 


土地 覆 被 类 别 样本 数量 /个 栅 格 数量 /个 
*I 100 12937 
灌木 102 9180 
草地 132 11223 
ARK 115 8511 
针 叶 林 157 20567 
Duis 109 16459 
针 阔 混交 林 48 5769 
建设 用 地 76 8360 
裸 地 103 8652 
水 域 157 15563 


似 表 示 地 表 波 形 信 号 。 对 其 求 差 得 到 的 波形 幅度 
(Loa ) 可 以 表示 GLAS 提 取 冠 层 高 度 。 
L. = d_gpCutRngOff - d. SigBegOff (1) 

在 平坦 地 区 , L, 可 以 近似 表示 实际 的 森林 冠 
层 高 度 ,但 在 山区 ,由 于 地 形 起 伏 影 响 ,GLAS 波形 
会 发 生变 化 ,由 此 引起 的 偏差 值 Bias 会 直接 影响 森 
林 冠 层 高 度 的 估计 结果 (图 4)。 

本 文采 用 1 种 基于 几何 模型 的 地 形 校 正方 法 ”， 
对 Bias 进行 量化 : 
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d d 


注 :及 :为 实际 冠 层 高 度 (m);0 为 地 表 坡 度 (°);d 为 GLAS 脚 印 直径 (m)。 
图 4 不 同 地 形 条 件 对 GLAS 提取 和 森林 冠 层 高 度 的 影响 
Fig.4 Effects of different topographical conditions 


on GLAS-derived canopy height 


L.a = H+ Bias (2) 


Bias =d x tan 9+ FWHM (3) 


式 中 : 及 .表示 实际 冠 层 高 度 (m) ; Bias 为 偏差 值 ;d 
表示 GLAS 足 迹 直 径 , 本 研究 中 取 值 为 55 m; ORR 
足迹 所 在 的 坡度 等 级 的 平均 坡度 值 (*);c 表示 光速 
[m * (ns) ;FWHM(Full width at half maximum ) 表 示 
激光 脉冲 宽度 ,为 6 nso 

由 于 GLAS 数 据 与 样 地 调查 数据 以 及 森林 资源 
清查 数据 间 存 在 约 10 a 的 时 间 差 ,因此 ,我 们 引入 平 
移 系数 a 和 5 对 此 模型 进行 改进 ,并 采用 通用 全 局 优 
化 算法 (Levenberg marquard,LM) 基 于 GLAS 脚 印 点 
实地 调查 数据 拟 合 a 和 2 的 最 优 估计 值 ,从 而 使 
GLAS 提取 的 冠 层 高 度 更 接近 研究 区 实际 的 冠 层 高 
度 。 模 型 改进 后 从 GLAS 数 据 中 提取 的 冠 层 高 度 为 : 

Haas =La lax (dxtand+ CI 4 (4) 


2.2.3 GLAS 脚 印 点 碳 密度 估算 ”人 研究 区 典型 林 分 
调查 共 建 立 针 叶 林 样 地 42 个 , 阔 叶 林 样 地 38 个 , 针 
阔 混 交 林 样 地 18 个 。 按 照 森 林 类 型 统计 各 样 地 的 
平均 树 高 (7) 和 平均 胸径 (D)。 借 鉴 地 区 物种 的 异 
速生 长 方程 ( 表 2) ,随机 选取 每 种 类 型 森林 70% 的 
样 地 调查 数据 , 拟 合 HD、DH IgH 1gD FillgD°H iX 6 
个 自 变量 与 地 上 生物 量 (AGB ) 的 函数 。 同 时 以 平均 
Ws CD) 2 ELE Rt ,以 平均 树 高 (万 ) 为 自 变 量 , 拟 合 
D lgD ,D? 5j H gH PARAR., li FELIS 
30% 的 样 地 数据 交叉 验证 各 函数 的 拟 合 效果 ,以 模 
型 决定 系数 尼 判 定 拟 合 效果 的 优 劣 。 

因 利 用 与 树种 相关 的 碳 含量 转换 系数 可 以 将 
地 上 生物 量 转换 为 地 上 碳 储量 时 得 到 最 准确 的 佑 
计 ' 引 。 因 此 ,对 进行 实地 调查 的 97 个 GLAS 脚 印 点 ， 
根据 表 2 选 择 特定 的 碳 含量 转换 系数 计算 其 碳 密 


度 。 对 于 剩余 的 176 个 GLAS 脚印 点 ,根据 区 域 相 关 
WER ,分 别 将 针 叶 林 、 闪 叶 林 和 针 冰 混交 林 的 地 
上 生物 量 按照 0.4927 .0.4682 和 0.4682 的 碳 含量 
换 系 数 转化 为 地 上 碳 储量 ,并 结合 森林 资源 清查 数 
据 中 记录 的 林 分 密度 计算 脚印 点 地 上 碳 密度 。 
2.2.4 研究 区 森林 地 上 碳 密度 制图 由 于 GLAS 脚 
印 点 的 不 连续 性 ,需要 与 辅助 数据 结合 以 获得 空间 
连续 的 森林 地 上 碳 密度 分 布 。 将 激光 脚印 点 的 碳 
密度 外 推 到 整个 研究 区 域 的 方法 通常 分 为 两 类 : 参 
数 化 算法 和 非 参数 算法 。 参 数 化 算法 ,一 般 使 用 具 
有 空间 连续 性 的 光学 遥感 数据 建立 基于 像素 的 生 
物 量 / 碳 储量 与 光学 指数 间 的 回归 模型 " ,但 实际 
上 ,森林 生物 量 受 多 种 因素 的 影响 ,其 与 遥感 数据 
的 关系 很 难 用 简单 的 线性 或 非 线 性 模型 表达 '。 
此 , 越 来 越 多 的 研究 选择 使 用 非 参 数 化 算法 将 激光 雷 
达 脚 印 点 生物 量 / 矶 储量 外 推 到 整个 研究 区 域 *7。 

本 人 研究 采用 非 参 数 化 算法 最 大 炉 模 型 (Max- 
Ent) 获 取 人 研究 区 森林 地 上 碳 密度 分 布 图 。 最 大 信 模 
型 是 基于 最 大 炉 原 理 用 于 预测 物种 潜在 地 理 分 布 
的 预测 模型 中 。 运 行 时 首先 确定 物种 己 知 点 的 分 
布 区 域 ,寻找 限制 物种 分 布 的 约束 条 件 即 环境 变 
i ,构建 约束 集合 ,根据 物种 实际 分 布点 和 环境 变 
量 得 出 预测 模型 ,再 对 模拟 目标 物种 在 目标 区 域 的 
可 能 分 布 进行 预测 ,以 所 有 满足 已 知 约束 条 件 集合 
HE od A AEE o 

运行 模型 时 以 研究 区 环境 及 地 形 数据 共计 23 
个 图 层 (19 个 生物 气候 变量 图 .土壤 质地 空间 分 布 
数据 以 及 海拔 .坡度 和 坡 向 图 ) 作 为 环境 变量 , 以 
GLAS 脚印 点 和 样 地 调查 点 的 地 上 碳 密 度 作 为 训练 
数据 。 将 样 地 调查 点 和 GLAS 脚印 点 地 上 碳 密 度 
估算 结果 划分 为 12 个 等 级 区 间 :5~10 thm? , 10~ 
20 t- hm™?, 20~30 t+ hm, +++, 100~110 t+ hm, 110~ 
119 thm > ,分 别 将 每 一 等 级 碳 密度 的 样 地 调查 和 
GLAS 脚印 点 的 地 理 位 置 输入 以 构建 该 等 级 地 上 碳 
密度 分 布 的 约束 环境 集合 ,并 生成 该 等 级 地 上 碳 密 
度 在 每 个 栅 格 概率 值 分 布 ,并 使 用 以 下 函数 将 其 转 
换 为 每 个 像素 的 地 上 碳 密度 : 


n 
Yr 
a tel 
= n 
P" 


i=l í 
RP: AGC, EARE LIRE n d EE P 
数 ; P? BRAK n UGE PRU BA A 4 8 


AGC ity (5 ) 
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表 2 祁连山 地 区 相关 树种 异 速 方程 和 碳 含量 转换 系数 
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Tab.2 Species-specific allometric equations and carbon content conversion factor in Qilian Mountains 


树种 异 速生 长 方程 参考 文献 碳 含 量 转 换 系数 
青海 云 杉 W.-0.0478(D'H) ^ 王 金 叶 等 ” 0.5200" 
Picea crassifolia W.-0 0122(D H) 
Hoe 
W,- 0.2650(D^ H) " 
W=Wo+ Ws + W, 
祁 连 圆 柏 W.-02738(D'H) - 王 金 叶 等 7 0.4650"! 
Juniperus przewalskii W.=0.0061 (pay 
Br Ms 
W, =0.0042(D°H) 
W=W,+W,+W, 
侧 柏 W.-0.0427(D'H) ^ KEES" 0.5034" 
Platycladus orientalis W =0.0128 (pH) 
» 70. 
W,-00258(D'H) . 
W=W,;+W,+ W, 
油 松 Wen 程 堂 仁 等 " 0.5108" 
Pinus tabulaeformis W = em (pH) 
= 
W =e’ "pg 
L 
W= Ws + Ws + W 
a W.=0.0493(D"H) KEKE 0.5004" 
Quercus mongolica W.-0 0026(D°H) 
Bo 
W,=0.0119(D*H) 
W=W,+W,+W, 
山 杨 W.- ga D Hy" Fee (Go 0.4956"! 
Populus davidiana W = e500 (pen) 
B EE 
W= (pH) 
L = 
W=Wo+ Ws + W, 
Si FE W. = grey 程 堂 仁 等 6 0.5025"! 
Betula platyphylla W =e (poa) 
B =e. 
W= (pH) 
L x 
W= Ws + Ws + W, 


地 上 碳 储 量 密度 分 布 的 概率 ; C, 是 每 个 地 上 碳 密 度 
等 级 区 间 的 中 值 。 经 过 多 次 迭代 和 与 测试 数据 的 
交叉 验证 ,发 现 当 n=3 时 ,地 上 碳 密度 分 布 收敛 到 最 
优 值 。 男 外 ,使 用 每 个 碳 密 度 等 级 的 最 大 值 和 最 小 
值 来 创建 概率 分 布 图 以 包含 估计 的 上 下 界 。 最 后 
根据 研究 区 土地 覆 被 分 类 的 森林 范围 对 所 得 结 
进行 担 膜 提取 ,以 避免 保留 毁 林 和 灌 从 地 区 地 上 碳 
密度 值 而 对 研究 区 森林 地 上 碳 密 度 造成 高 佑 。 


3 结果 与 分 析 


3.1 GLAS 脚 印 点 森林 冠 层 高 度 估算 结果 
通过 运行 LM 算法 获得 各 坡度 等 级 平移 系数 a 


注 :两 为 干 生 物 量 (kg);W 为 村 3 生物 量 (kg); 极为 叶 生 物 量 (kg); 丈 为 地 上 生物 量 (kg ); 刀 为 胸径 (cm ) ;万 为 树 高 Am)。 


和 2 的 最 佳 估计 值 。 按 照 不 同 坡度 等 级 对 GLAS 脚 
印 点 提取 的 森林 冠 层 高 度 进行 校正 ,并 随机 选择 51 
个 GLAS 脚印 点 以 森林 资源 清查 数据 独立 验证 
GLAS 数 据 提取 的 冠 层 高 度 精度 及 地 形 校正 效果 ， 
地 形 校正 效果 如 图 5 所 示 。 

结果 表明 ,与 校正 前 相 比 ,通过 改进 后 的 地 
校正 模型 提取 的 冠 层 高 度 更 接近 1:1 拟 合 线 , 即 更 
接近 森林 资源 清查 数据 中 记录 的 冠 层 高 度 。 表 3 显 
示 了 各 坡度 等 级 地 形 校 正 模型 的 统计 精度 。 

由 表 3 可 知 , 地 形 校 正 后 ,各 坡度 等 级 地 形 校 正 
模型 的 让 可 以 解释 65%~67% 的 冠 层 高 度 的 变化 程 
度 ,与 校正 前 的 尼 相 比 变化 并 不 明显 ,表明 改进 后 


SN 
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— 1:18 s 校正 前 © 校正 后 
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30r (b) 5°~15° 坡 度 级 
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5 10 15 20 25 30 
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30r (d) 25°%~35° 坡 度 级 
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30[ (e) 35°~58° 坡 度 级 


GLAS 提 取 冠 层 高 度 /m 


10 
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5 地 形 校 正 前 后 不 同 坡度 等 级 GLAS 冠 层 高 度 估 计 精 度 对 比 


Fig. 5 Comparison of GLAS derived canopy height before and after topographic correction in different slope gradient classes 


表 3 不 同 坡度 等 级 地 形 校正 前 后 GLAS 提取 
冠 层 高 度 的 精度 比较 
Tab.3 Accuracy comparison of GLAS derived 
canopy heights before and after topographic 
correction in different slope gradient classes 


坡度 等 级 R RMSE/m 
校正 前 校正 后 校正 前 校正 后 
0°~5° 0.62 0.66 2.83 1.84 
59-15? 0.64 0.65 3.32 2.23 
159-25? 0.66 0.67 6.74 3.51 
259-35? 0.63 0.65 8.00 3.91 
35°~58° 0.65 0.67 12.59 3.48 


注 : 尼 为 模型 决定 系数 ;RMSE 为 模型 均 方 根 误差 。 


的 地 形 校 正 模型 并 没有 改变 数据 的 拟 合 优 度 。 男 
外 ,地 形 校正 前 , 随 着 坡度 的 增加 , La A9 RMSE 从 
2.83 mm 增加 到 12.59 m。 地 形 校正 后 ,不 同 坡度 等 级 
的 Hos 的 RMSE 稳 定 在 1.84~3.91 m 之 间 , 因 此 ,我 
们 认为 改进 后 的 地 形 修正 模型 性 能 良好 , 且 对 较 高 
坡度 等 级 的 地 形 修正 效果 更 加 明显 。 
3.2 GLAS 脚印 点 地 上 碳 储量 密度 估算 结果 
根据 研究 区 典型 林 分 调查 结果 和 树种 异 速生 
长 方程 ,建立 不 同类 型 森林 的 地 上 生物 量 估算 模 
型 。 拟 合 过 程 中 发 现 , 仅 加 入 树 高 信息 时 模型 拟 合 
精度 最 低 ,而 同时 加 入 树 高 胸径 信息 时 模型 估 测 精 
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度 最 高 。 最 终 针 叶 林 和 阔 叶 林 的 地 上 生物 量 估算 
模型 以 虞 函数 拟 合 效 果 最 佳 ,其 模型 决定 系数 尼 分 
别 为 0.979 和 0.983 , 针 疗 混 交 林 的 地 上 生物 量 估 算 
模型 以 二 次 项 困 数 拟 合 效果 最 佳 , 尼 为 0.987。3 
种 不 同 森 林 类 型 的 地 上 生物 量 估算 模型 如 式 (6)~ 
式 (8): 


AGB( 针 叶 林 )=0.127x(D DD (6) 

AGB (RAF HK) = 1.444 x (D? H^? (7) 

AGB (#t AYE Ze #4) = 275.42 x (Ig D HY — 
1808.383 x lg D^ H + 3078.138 

由 于 未 进行 实地 调查 的 GLAS 脚印 点 缺少 树木 

胸径 数据 ,另外 按照 不 同 森 林 类 型 对 其 树 高 和 胸径 

的 关系 进行 拟 合 。 最 终 针 叶 林 树 高 与 胸径 间 的 关 

系 方程 以 $ 函数 拟 合 效果 最 佳 , 尼 为 0.787; 阀 叶 林 

以 二 次 项 函数 拟 合 效果 最 佳 ,RR 为 0.778; 针 阔 混 交 

林 以 复合 函数 拟 合 效果 最 佳 , 尼 为 0.734。3 种 不 同 

类 型 森林 树 高 胸径 关系 如 式 (9)~ 式 (11): 

lg D (EFM Bh) =e tmm (9) 

D*( Kalk) = 0.0034 x (^Y + 1.127 x H? + 384.96 (10) 

D? (EF RAYEZS PR) = 77.71 x 1.175" (11) 

利用 典型 林 分 样 地 调查 数据 交叉 验证 以 上 模 

型 单 木 水 平 的 估计 精度 : 针 叶 林 、 阔 叶 林 和 人 针 痊 混 


(8) 


交 林 森林 地 上 生物 量 估算 的 均 方 根 误差 RMSE 分别 
为 17.88 kg, 7.13 kg 和 4.65 ke, aE TH TAY RMSE 分 
HIJA 5.2 cm .6.98 cm 和 4.04 cmo 

按照 以 上 生物 量 估 算 模 型 及 树 高 胸径 关系 计 
算 位 于 不 同 森 林 类 型 CLAS 脚印 点 内 的 地 上 生物 
量 , 并 根据 各 森林 类 型 的 碳 含 量 转换 系数 以 及 森林 
调查 数据 中 记录 的 林 分 密度 得 到 森林 区 域 GLAS 肢 
印 点 森林 地 上 碳 密 度 , 估 算 结 果 如 图 6 所 示 。 
3.3 研究 区 森林 地 上 碳 密度 分 布 估算 结果 

通过 MaxEnt 模 型 将 GLAS 脚印 点 森林 地 上 碳 密 
度 外 推 至 整个 研究 区 域 ,得 到 2018 年 祁连山 国家 公 
30 m 分 辩 率 的 森林 地 上 碳 密度 空间 分 布 如 网 7 
所 示 。 

为 验证 估算 结果 的 准确 性 ,从 与 GLAS 脚印 点 
及 样 地 调查 点 至 少 相 距 1 km 的 森林 资源 清查 数据 
中 随机 选取 181 个 样本 点 独立 验证 估算 精度 ,得 到 
本 次 森林 地 上 碳 密 度 估算 均 方 根 误差 RMSE 为 
18.946 thm?。 此 结果 满足 REDD+ 计 划 MRV(Mea- 
suring, reporting, and verification ) 指南 的 要 求 , 即 利 
用 卫星 遥感 数据 估算 森林 碳 储量 时 ,应 将 佑 算 值 与 
实地 测量 值 之 间 的 误差 控制 在 20 thm ?或 20% 
以 内 。 

结果 显示 ,2018 年 祁连山 国家 公园 森林 地 上 碳 
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图 6 森林 区 域内 GLAS 脚印 点 森林 地 上 碳 密 度 估算 结果 


Fig.6 Estimation of forest aboveground carbon density of GLAS footprints in forest area 
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图 7 2018 年 祁连山 国家 公园 森林 地 上 碳 密度 分 布 图 
Fig. 7 Forest aboveground carbon density map of Qilian Mountains National Park in 2018 


密度 主要 分 布 在 20.67~62.62 thm ”之 间 , 最 高 值 为 
119.18 t- hm^, 平均 森林 地 上 碳 密 度 为 40.72+6.72 
thm ,和 森林 地 上 总 碳 储量 为 28.58+4.72 Tg(1 Tg = 
10 kg)。 其 中 甘肃 省 境内 森林 地 上 碳 储量 为 23.67 
Tg, 远 高 于 青海 省 (5.21 Tg) ,甘肃 省 森林 地 上 碳 密 
度 平 均值 为 44.90 t- hm? ,也 高 于 青海 省 地 上 碳 密度 
平均 值 28.71t:hm?*。 不 同类 型 森林 中 ,和 针 叶 林 的 平 
均 地 上 碳 密 度 最 高 (50.69 thm?) , 接 下 来 依次 是 阅 
叶 林 (41.42 t* hm?) FFP SHE SAK (38.06 t* hm?) , 针 
叶 林 BRD RAR BE TEC 36 AS EBK tae 23 AR 
连 山 国家 公园 森林 地 上 总 碳 储量 的 53.28% 16.57% 
Fil 30.15%. 

由 于 植被 分 布 主要 受气 候 变 化 的 驱动 特别 是 
温度 和 降水 的 影响 ,而 山区 高 原 复杂 的 地 形 对 气候 
因子 的 空间 格局 作用 较 大 ,因此 ,我 们 通过 地 学 
统计 进一步 对 位 于 不 同 海拔 等 级 (1770~2770 m、 
2770~3770 m,3770-4770 m ,4770~5740 m) .坡度 等 
2K (09-59 5°~ 15°, 15°~25° ,25°~35° All 35° ~58° ) WK 
坡 向 ( 阳 坡 、 阴 坡 、 半 阳 坡 以 及 半 阴 坡 ) 的 森林 地 上 
碳 储量 进行 分 析 。 结 果 显 示 ,对 不同 的 海拔 等 级 而 
Ti ,海拔 2770~3770 m 地 区 森林 地 上 碳 储量 最 大 , 约 
占 研 究 区 森林 地 上 总 碳 储量 的 64.47% ,其 次 是 海拔 


1770-2770 m 地 区 (21.05% ) 以 及 海拔 3770~4770 m 
地 区 (13.99% ) ,海拔 4770~5740 m 地 区 基本 没有 森 
林地 上 碳 储量 分 布 。 对 不 同 的 坡度 等 级 而 言 BR 
35°~58° 地 区 森林 地 上 碳 储 量 较 少 外 ,其 余 坡度 的 森 
林地 上 碳 储 量 基本 相同 。 对 不 同 坡 向 而 言 , 阴 坡 及 
半 阴 坡 的 森林 地 上 碳 储量 基本 为 阳 坡 及 半 阳 坡 的 
2 倍 。 


4 讨论 与 结论 


4.1 讨论 

本 研究 结果 小 于 Wagner 等 ”通过 采样 法 获得 
的 祁连山 区 东北 缘 青 海 云 杉林 地 上 碳 密度 (82.70 
t* hm?) ,这 可 能 是 因为 该 区 域 树 木 生长 状况 较 好 且 
青海 云 杉林 为 区 域 优 势 树种 。 本 研究 结果 与 Zhao 
等 6 结合 LPJ-GCUESS 和 HASM 模 型 对 中 国 西北 地 区 
森林 地 上 碳 密 度 的 估计 结果 (36.10~64.17 thm?) Be 
近 。 将 本 研究 结果 与 其 他 山区 和 森林 地 上 碳 密度 佑 
算 结果 进行 比较 ,本 研究 祁连山 区 森林 地 上 和 平均 碳 
密度 估算 结果 (40.72+6.72 t* hm?) SJA XEF A 
然 保 护 区 "森林 地 上 平均 碳 密度 (53.14 t hm?) 以 
及 内 蒙古 大 兴安 岭 生 态 站 “森林 地 上 平均 碳 密 度 
(38.29 thm 7) 接 近 , 说 明 祁 连 山 与 其 他 重要 的 山地 
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自然 保护 区 一 致 ,是 重要 的 高 海拔 森林 固 碳 场 所 。 

本 研究 依据 GLAS 数据 提供 的 森林 垂直 结构 信 
息 估 测 森 林地 上 碳 密度 。 目 前 从 GLAS 数 据 中 提取 
森林 和 冠 层 高 度 的 方法 主要 分 为 两 类 :第 一 类 是 将 第 
一 个 波峰 和 地 面 波峰 之 间 的 距离 作为 植被 冠 层 的 
高 度 值 ,而 地 面 波峰 通常 被 认为 是 最 低 处 的 高 斯 峰 
值 或 最 低 2 个 高 斯 峰值 中 较 强 的 高 斯 峰值 汪 ;第 二 
类 是 根据 波形 范围 或 其 他 波形 参数 通过 统计 分 析 
拟 合 植被 高 度 ，。 比 较 而 言 ,第 一 类 方法 较为 简 
单 ,可 以 灵活 地 应 用 于 不 同 的 地 区 ,特别 是 结合 
DEM 数据 等 其 他 辅助 数据 时 ,能 够 较为 直观 的 消除 
地 形 对 GLAS 波 形 的 影响 ,同时 该 方法 也 被 诸多 研 
究 证 实在 山区 森林 冠 层 高 度 估算 方面 具有 较 好 的 
性 能 了。 因此 ,本 研究 采用 第 一 类 方法 ,结合 改进 
后 的 地 形 校 正 模 型 从 GLAS 数 据 中 提取 和 森林 冠 层 高 
度 。 另 外 ,本 研究 选择 已 经 过 统一 滤波 和 拟 合 处 理 
的 GLA14 产品 而 非 GLAOL 产品 ,可 以 减少 处 理 
GLA01 产品 中 记录 的 原始 波形 时 的 一 些 侦 然 误差 。 

本 研究 采用 MaxEnt 模 型 将 GLAS 脚印 点 森林 地 
上 碳 密 度 外 推 至 整个 研究 区 域 。MaxEnt 是 1 种 基于 
概率 框架 从 一 组 不 完全 信息 中 生成 预测 或 推断 的 
通用 算法 中 。 它 依赖 于 不 完全 经 验 概率 分 布 的 假 
设 , 在 一 定 的 环境 约束 下 ,该 假设 可 以 近似 为 具有 
最 大 信 的 概率 分 布 (MaxEnt 分 布 ), 并 且 该 分 布 接近 
于 潜在 的 地 理 分 布 ”。 考 虑 到 研究 区 内 物种 分 布 
主要 受气 候 因 素 的 控制 且 复杂 地 形 对 物种 分 布 也 
有 很 大 影响 ,本 研究 采用 MaxEnt 模 型 基于 GLAS 脚 
印 点 和 样 地 调查 点 的 地 上 碳 密 度 及 其 分 布 的 气候 、 
土壤 及 地 形 数据 作为 训练 数据 ,模拟 接近 实际 情况 
的 森林 地 上 碳 密度 分 布 。 

已 有 人 研究 证 实 利 用 遥感 数据 估算 森林 地 上 生 
物 量 / 碳 储量 的 主要 误差 源 是 不 同 数据 源 之 间 的 空 
间 异 质 性 所 。 而 本 研究 的 空间 异 质 性 问题 主要 涉 
及 GLAS 足 迹 的 位 置 精 度 以 及 所 有 数据 空间 分 辩 率 
的 一 致 性 。 因 此 ,我 们 首先 依据 GLA14 产 品 中 记录 
的 高 程 用 DEM 数据 检查 GLAS 足迹 定位 精度 。 其 
次 ,我 们 选择 仅 来 自 3 号 激光 器 的 GLAS 数 据 并 设计 
了 与 其 分 辩 率 匹配 的 样 地 调查 方案 ,同时 还 将 所 有 
影像 数据 分 辨 率 重 采样 到 30 m, 以 最 大 程度 的 消除 
由 空间 异 质 性 引起 的 误差 。 

本 研究 的 不 确定 性 包括 以 下 几 个 方面 :(1) FE 
地 调查 时 的 测量 误差 ;(2) EH BE BEATE AL A AY 


类 ;(3) 冬季 获取 的 GLAS 数 据 导致 的 由 落叶 树种 或 
冠 层 积 雪 引 起 的 森林 冠 层 高 度 估 计 误 差 。 
4.2 结论 

本 人 研究 基于 CLAS 数 据 、Landsat 0LI 数 据 以 及 
样 地 调查 数据 ,在 建立 不 同类 型 森林 地 上 生物 量 佑 
算 模型 的 基础 上 ,利用 GLAS 数 据 提取 的 森林 冠 层 
高 度 对 其 脚印 点 森林 地 上 碳 密度 进行 估算 ,并 通过 
MaxEnt 非 参 数 化 方法 将 脚印 点 森林 地 上 碳 密 度 外 
推 到 整个 研究 区 域 。 结 果 表 明 : 

(1) 从 GLAS 数 据 中 提取 和 森林 冠 层 高 度 时 ,通过 
改进 后 的 地 形 校正 模型 提升 GLAS 数据 提取 的 精 
RE ,这 不 仅 为 GLAS 脚印 点 森林 地 上 碳 密度 的 准确 
估算 提供 了 保证 ,也 使 得 更 多 的 GLAS 数据 可 用 于 
后 续 的 研究 。 

(2) 所 得 的 祁连山 国家 公园 森林 地 上 碳 密度 佑 
算 结 果 满 足 REDD+ 计 划 MRYV 指南 对 利用 卫星 遥感 
数据 估算 森林 碳 储量 的 误差 要 求 ,可 为 监测 区 域 乃 
至 国家 尺度 的 森林 碳 储 量变 化 以 及 制定 可 持续 的 
森林 管理 措施 提供 依据 。 此 外 ,本 文 所 采用 的 方法 
在 山区 森林 碳 储 量 估算 方面 也 具有 较 大 的 潜力 。 
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Estimation of forest aboveground carbon density in Qilian 
Mountains National Park based on remote sensing 
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(1. College of Natural Resources and Environment, Gansu Agricultural University, Lanzhou 730070, Gansu, China; 


2. Land Use Research Institute, Gansu Agricultural University, Lanzhou 730070, Gansu, China) 


Abstract: An accurate, up-to-date, and spatially explicit estimate of the forest aboveground carbon (AGC) density 
is indispensable for understanding the global carbon cycle and implementing measures to mitigate climate change 
such as reducing emissions from deforestation and forest degradation (REDD+). Geoscience Laser Altimeter System 
(GLAS) data have proven to be a powerful candidate for estimating the forest AGC density distribution over large- 
scale areas because of its global-scale sampling strategy and freely available data. However, the influence of terrain 
slopes on the accuracy of the GLAS- derived canopy height often limits its utility in mountainous regions. 
Nevertheless, mountain forests play a key role in carbon storage because of the relatively small influence of 
human activities, and the carbon sequestration capacities of forests differ according to the region. In this study, 
we integrated GLAS data, Landsat 8 OLI imagery, and field inventory data to estimate the forest AGC density 
distribution of Qilian Mountains National Park in northwestern China, which is an essential ecological barrier and 
carbon sequestration site in western China. We aimed to address the underestimation of carbon storage by 
mountain forests in previous studies by improving the estimation accuracy of the GLAS-derived canopy height 
and in turn the forest AGC density in the study area. We first improved upon a physical terrain correction model 
to derive the canopy height from GLAS data according to five different slope gradient classes. The canopy height 
could be reliably estimated with a root- mean- squared error (RMSE) of 1.84-3.91 m on the basis of independent 
validation with forest inventory data. We subsequently calculated the AGC density of GLAS footprints located in 
different forest types as accurately as possible by using an aboveground biomass estimation model of different 
forest types and a carbon content conversion factor, as well as the stand density recorded in the forest inventory 
data. We took a nonparametric approach and applied a MaxEnt model to extrapolate AGC data at the GLAS 
footprint scale to the whole forest as extracted from Landsat 8 OLI imagery. We also made great efforts to 
spatiotemporally match all data sources in this study. The results showed that the average forest AGC density in 
Qilian Mountains National Park was 40.72-56.72 t- hm ? in 2018, and the total aboveground carbon storage was 
28.58+4.72 Tg. Forest vegetation in areas at 2770-3770 m above sea level had the largest carbon storage, and 
shaded slopes had significantly greater carbon storage than sunny slopes. The accuracy of the estimation results 
was independently verified by comparison with the forest resource inventory data and resulted in an RMSE of 
18.946 t* hm ^, which meets the measuring, reporting, and verification guidelines of REDD+. The results of this 
study can provide a basis for monitoring regional- and even national-scale changes in forest carbon reserves and 
for formulating sustainable forest management policies. Additionally, the method used in this study can 
potentially be applied to estimating the carbon storage of forests in mountainous areas. 

Key words: satellite LIDAR system; Landsat OLI; forest aboveground carbon density; non-parametric algo- 
rithm; Qilian Mountains 


